Anisotropic conductivity in superconducting NCCO 
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The low temperature behaviour of the in-plane and c-axis conductivity of electron-doped cuprates 
like NCCO is examined; it is shown to be consistent with an isotropic quasiparticle scattering rate 
and an anisotropic interlayer hopping parameter which is non-zero for planar momenta along the 
direction of the d x 2_ y 2 order parameter nodes. Based on these hypotheses we find that both, the in- 
plane and the c-axis conductivity, vary linearly with temperature, in agreement with experimental 
data at millimiter-wave frequencies. 

PACS numbers: 74.72.-h, 74.72. Jt, 72.10.-d 
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The pairing symmetry of electron-doped cuprate ma- 
terials has been the subject of renewed interest in the 
last year. Previously existing experimental data were 
interpreted to be consistent with an s-wave order pa- 
rameter 0. In particular penetration depth measure- 
ments M, Raman scattering 0, and tunneling data Q 
were explained in this way. Surprisingly, however, recent 
microwave experiments and phase-sensitive tricrys- 
tal experiments ||] have provided evidence in favor of a 
d 2 ,2_j ( 2-wave pairing symmetry. 

In this report we examine the low temperature be- 
haviour of the anisotropic conductivity in superconduct- 
ing electron doped cuprates. Using a tilted film geom- 
etry, Pimenov et al. |12f have recently measured simul- 
taneously <7 C and a a b (c-axis and in-plane conductivity) 
on Nd2- :c Ce :I ;Cu04 (NCCO) samples. By measuring the 
transmission through thin films as well as the transmis- 
sion phase shift with an interferometer, both complex 
conductivities were determined in the submillimiter fre- 
quency range. Their results show a characteristic temper- 
ature dependence of a c and a a b in the superconducting 
state, varying both linearly with T almost up to T c . Fig. 
|l| shows an example of these experimental data |l2] for a 
NCCO film with T c « 20 K. Fits to the low temperature 
data are added to underline the linear T dependence of 
both conductivities. In what follows we show how this 
linear T dependence is consistent with a d x 2_ a 2-wave or- 
der parameter under certain assumptions for the scatter- 
ing rate (T^) and the c-axis hopping amplitude (tj_(k)). 

For tetragonal cuprate compounds the interlayer c-axis 
hopping amplitude has the form H 



t±(k) = — 1± cos(k z d)(cos(k x a) — cos(k y a)) i 



(1) 



where k x , k y are the components of the in-plane momen- 
tum, k z the momentum along the c axis and a and d are 
the corresponding lattice constants. From now on a will 
serve as a length unit and is set to 1. The form of tj_(k) 
in Eq. (Q) results from the hybridization between the 
bonding oxygen 2p and copper 4s orbitals in each Cu02 
plane which gives rise to the nodal structure of Eq. (Q) . 
A small deviation from the tetragonal structure leads to 
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FIG. 1. Experimental data for the real part of the c-axis 
(cr c ) and the in-plane conductivity (<j a b) of a NCCO film with 
T c = 20 K at the frequency of v = 8 cm -1 (taken from Ref. 
[11]). The solid lines are linear fits to the low temperature 
part of the data. 



a finite value of t± (k) also along the planar nodal direc- 
tions. Similarly, isotropic impurity scattering is expected 
to effectively contribute to the c-axis transport. In this 
report we consider therefore the following extended form 
fortx(k) 

tj_(k) = — t±[cos(k z d)(cos(k x ) — cos(k y )) 2 + c]. (2) 

We are interested in the real part of the dc limit of the 
longitudinal conductivity both for the in-plane and the 
c-axis case. The conductivity, labeled by er^, where [i is 
either zon, is given by: 

Im D^(m + iO+) 
a = Re cr w = lim ^ (3) 

where D^^iQ) is the current-current correlation function 



lim 



e inr (T T j M (q,r)i M (-q,0))dr (4) 



with j M (q, t) the paramagnetic current density operator. 
Note that the usual notation a c corresponds to a zz and 



1 



Cab to a xx . In Eq. (Q) /? = 1/T, T the temperature, T T 
the imaginary time ordering operator and fl is a bosonic 
Matsubara frequency. For the purpose of comparing with 
experimental data in the submillimcter frequency range 
presented in Fig. ^ the zero frequency limit in Eq. (j^) is 
appropriate. 

Neglecting vertex corrections we rewrite D PAl (if2) as 



(•5) 



ja;( z )(k) = e/h(dt;k/dk x ( z j) is the current vertex function 
and G(k, iui) is the renormalizcd BCS Green's function 
in Nambu space 



G(k,iw) 



(iU))T + CkT 3 + A k Ti 



£ 2 



(0) 



where cD = uj — Sq, £k = £k + S3 and Ak = Ak + Si. 
£k is the band dispersion measured with respect to the 
chemical potential, £k = £k~A*, Tj are the Pauli matrices, 
and Ak is the superconducting gap that we take in a 



,,2 -wave form 



A k = ^(cos(fe x ) - cos(fcj,)). 



(7) 



Si are the components of the self energy matrix S = 
£i=o ^i T i assumed to dominantly arise from impurity 
scattering in the low frequency and low temperature limit 
that we consider here. Using the spectral representation 
for G and summing over the internal Matsubara frequen- 
cies the real part of the dc conductivity results as 



^ = TTEtfM / &M^Tr{ImG(k, £ ) ; 



(8) 



Im G(k, e)} 



where f(x) — (exp(ai) + 1) 1 is the Fermi function. In 
the following we assume that the band dispersion £k de- 
pends only on k z and k = 
momentum sum as 



k 2 and we evaluate the 



1^ - at(o) r dk z r deb f d£ k . 

1 J — TT J <J 



(9) 



cb is the azimuthal angle around the cylindrical Fermi 
surface (FS) and the density of states is approximated by 
its value at the FS N(£ k ) « N(0). All the functions in 
Eq. (|8|) are then parametrized along the FS as functions 
of cb, k z and £k- Note that in what follows we neglect 
the effect of Si, i.e. any gap renormalization, since we 
are interested in the effect of the nodal structure of a 
pure d 2 ,2_ ?/ 2-wave gap. For the purpose of calculating 
the dc conductivity limit it is furthermore justified to 



ignore Im S3 which is negligibly small in the very low 
frequency limit Jl(|. We set So = — iT where T is an 
isotropic impurity scattering rate on the FS. We note that 
an anisotropic form for So leads to higher order terms in 
the subsequent low temperature expansion |Tl| ] but leaves 
unaffected the leading low temperature behaviour. 

With a d-wave gap A^ = A cos(2</>) and t±(cb) — 
—t±[cos(k z d) cos 2 2cb + c] the result after integration over 
£k and k z is written in the form: 



■fifj, 



Re 



dc 



{c + iTf 



df(e/T) 
de 



(10) 



r[(e + *r) 2 - a 2 ]3/ 2 



Here we have defined a xx = 2he 2 N(0)v F with vf the 
Fermi velocity in the ab plane, a zz ■ 



2 N{0)t 2 ± d 2 /h, 



,(4>) = 1 and u zz (4>) = (c 



The energy integral in Eq. ( |10| ) can be divided in two 
parts, the first with |e| £ [0,T] and the second with |e| G 
[r, 00]. To analyze these two contributions separately we 
introduce the following notation 



(11) 



where, comparing Eq. (|Io| ) and Eq. (|TT|), 

f 2 * (e + iT) 3 - eA 2 

A w( £ ) = ~ a W J Q & ^W Re r[(e + , r)2 „ A2] 3/2 - 



(12) 



Eq. ([!(]) is thus rewritten as 



00 df(clT) 

de ^l/—LA^( e )e(-\e\+T)+ 



dc 



dc 

df(e/T) 
dc 



(13) 



A w (e)6(|e|-r) 



where Q(x) is the step function. Since we are interested 
in the low temperatures behaviour we perform a Som- 
merfeld like expansion of the first integral, leading to a 
dominant T independent contribution a' „„. The second 
integral is then evaluated using an expansion in e/T for 
the function A AlM (e) ||. 

The resulting leading terms for the longitudinal con- 
ductivity in powers of T/Ao are given by 

, , , e 2 4^(ln2)jV(0)^ 2 c 2 2 
°zz = a o,zz + T (J/A ) + G(J ) 



e 2 87r(ln 2)^(0)^71 



(T/A ) + Q(T 2 ). 



(14) 



These expansions are valid in the limit T < Ao expected 
to apply for NCCO [pi. 



2 



0, T 

Explicitly they are given by 



<7q jMA1 are the limits oj 



of the conductivities. 



4e 2 N(0)tld 2 c 2 
8e 2 N(0)v 2 F h 



'zz,0 



'xx,0 



(15) 



From Eq. ( jl5| ) we note that u' xx is independent of impu- 
rity controlled parameters such as T and possibly c. This 
fact is consistent with the work by Lee psj ], in which a 
universal limit of the optical conductivity for uo — > and 
T — > is found, i.e. independent of impurity scattering. 

The linear temperature dependence of <j' zz and a' xx in 
Eq. ( |l4| ) is clearly consistent with the experimental data 
on NCCO Q . In addition we may ask whether Eq. jl^ ) 
is also quantitatively consistent with the data. We have 
argued that c originates from small effects like deviations 
from the ideal tetragonal crystal structure or isotropic 
impurity scattering, and so we expect t±(<j) — 7r/4) <C 
t±(4> — 0), or equivalently c <C 1. We estimate c from 
the ratio r of the slopes of the two conductivities a' and 



i.e. 



' zzS) 



J xx,0 



2h 2 v 2 F 



(16) 



From the data on NCCO shown in Fig. [I] we find r » 
3 • 1CT 4 . Using d ~ 3 • 10~ 10 m, v F « 2.2 • 10 5 m/s |§ 
and £j_ « 10 meV the resulting estimate for c is 



10" 



(17) 



verifying our original assumptions. 

We point out that the result obtained in Eq. (|l4|) is 
due to the d x 2_ y 2--wave form for the superconducting gap 
combined with the modified c-axis hopping amplitude Eq 
(|f). A different symmetry for the order parameter leads 
to a very different low temperature behaviour for the 
conductivities. To this extent it is useful to compare 
these results with those calculated for an s-wave gap. 
Substituting = Ao = const, in Eq. ([Toj ) the low 
temperature expansion is now given by 



It is important to note that the only difference between 
our assumptions and the ones used for YBCO is the pres- 
ence of the constant c in t± (k) that allows us to obtain a 
low temperature behaviour of a compatible with NCCO 
data. 

In conclusion we have shown how, with an isotropic 
scattering rate and an intcrlayer hopping integral non- 
zero along the Brillouin zone diagonals, it is possible to 
obtain a linear T dependence for the low temperature be- 
haviour of the conductivities at millimeter wavelengths. 
This is indeed observed experimentally and supports a 
d 2 ,2_ l) 2-wave pairing symmetry for electron-doped mate- 
rials like NCCO. 

This work was partially supported by the Deutsche 
Forschungsgemeinschaft through SFB 484. 
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Xo / t (v / a u 7t + o(Vt7a^) 

A °/ T (VAo7r + 0(v^yAo")) 



(18) 

where c zz and c xx are two constants. This result is clearly 
incompatible with the experimental data. 

A similar analysis as presented above has been pre- 
viously applied to cuprates superconductors like YBCO 
fill, p[ with a different form for the interlayer hopping am- 
plitude, namely Eq. (Q) has been used with c = 0. This 
leads to a different low temperature behaviour for a c [[ll| , 
namely T 3 , that was shown to be compatible with experi- 
mental data on optimally doped YBCO and BSCCO O . 
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